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ABSTRACT 

^ , By decomposing the damping tail of CMB anisotropics into a series of transfer functions 

' representing individual physical effects, we provide ingredients that will aid in the reconstruction 

. of the cosmological model from small-scale CMB anisotropy data. We accurately calibrate the 

I model-independent effects of diffusion and reionization damping which provide potentially the 

■ most robust information on the background cosmology. Removing these effects, we uncover 

. model-dependent processes such as the acoustic peak modulation and gravitational enhancement 

I that can help distinguish between alternate models of structure formation and provide windows 

• into the evolution of fluctuations at various stages in their growth. 

o , 
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' 1. Introduction 

Much effort is being expended to measure the angular power spectrum of the cosmic microwave 
background (CMB) anisotropy on increasingly smaller angular scales. For many types of models for structure 
formation, the spectrum can be predicted to a precision of about 1% ( Hu et al. 1995| ), raising the hope that the 



cosmological parameters that are the input to these calculations can be extracted to comparable precision 



(see e.g. Jungman et al. 1996). The "inverse problem" of reconstructing the model given a spectrum is 
less well understood than the "forward problem" of predicting it given the model. For this purpose, it 
is important to assess the generation of anisotropics in a manner that is not tied to any given model for 
structure formation. From the theory of anisotropy formation, we know that CMB fluctuations suffered 
causal processing and damping of the primordial signal. In this paper, we numerically calibrate such effects, 



extending and improving upon prior work (Hu & Sugiyama 1995ab, 1996, hereafter HSa, HSb, HSc and Hu 



fc White 1996b) 



A particularly fruitful way to visualize the CMB spectrum, and one that provides a framework for the 
inverse problem, is as a product of transfer functions representing individual physical effects. The spectrum 
is then constructed out of physical elements rather than a model-dependent parameterization. Conceptually, 
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the evolution of CMB anisotropies processes primordial metric or gravitational potential perturbations into 
features observable in the spectrum today (see e.g. Bond 1996 , Hu fc White 1996b ). Since the evolution 
obeys linear perturbation theory, its effects are described by a series of transfer functions that multiply the 
underlying perturbations. The form of these functions depends on the cosmological model, not only for the 



background expansion and thermal history (see e.g. Bond et al. 1994, Seljak 1994, Hu, Sugiyama & Silk 



1996) but also for structure formation (see e.g. Crittenden & Turok 1995, Magueijo et al. 1996, Durrer et 



al. 1996, Hu, Spergel, & White 1996). By decomposing the evolution into functions representing separate 



physical effects, we can isolate portions of the anisotropy spectrum that are the most sensitive to particular 
aspects of the cosmological model. 



In particular, processes that damp CMB anisotropies, photon diffusion (Silk 1968) and rescattering 



(Efstathiou & Bond 1987), depend mainly on the background parameters and little on the perturbations 



that form structure in the universe. In §|^, we isolate these effects in a numerical treatment. From this 
damping calibration, we produce convenient fitting formulae that accurately describe the behavior of the 
diffusion and reionization damping transfer functions, or envelopes, directly in anisotropy multipole space. 
In we illustrate the reconstruction process by testing it with known models within the cold dark matter 
(CDM) scenario. By removing the model-independent effects of damping, one uncovers important model- 



dependent effects such as the baryon-drag modulation of the peaks (HSa), the potential envelope that 
describes gravitational driving of acoustic oscillations (HSc), and the regeneration of anisotropies during 



reionization (Sunyaev & Zel'dovich 197C, Kaiser 1984) 



In the context of currently popular models, recovery of these signatures will help distinguish between such 



possibilities as an inflationary or cosmological defect origin of fluctuations (Crittenden & Turok 199£ , Durrer 



et al. 1996 , Hu fc White 1996b ). The effects of damping are also intrinsically interesting because they provide 
the most model-independent probes of the background cosmology. We also consider how diffusion damping 
can be used to constrain the curvature of the universe and reionization damping to determine the redshift and 
extent of reionization in the universe. In this way, the study of effects in the damping tail of CMB anisotropies 
presented here will aid in the future reconstruction of the cosmological model from the anisotropy data. 



2. Damping Calculation 



Damping processes which affect CMB anisotropies provide the most mo del- independent information 
available in the spectrum and allow constraints on cosmological parameters such as the curvature and the 
thermal history of the universe. Furthermore, these universal effects obscure the model- dependent signatures 
that are useful to determine the mechanism for structure formation in the universe and the ultimate source 
of density perturbations. 

For both these reasons, an accurate calibration of damping effects is desirable. In this section, we begin 



with the formalism necessary to describe them (§2.1) and simple approximations to help understand their 
nature (§2.2). We then turn to numerical calibration of these effects (§2.3). Finally, we give convenient 



fitting formulae to their effects on the anisotropy power spectrum that encapsulate these results (§2.4) 
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2.1. Boltzmann Formalism 



In this section, we provide the formahsm for the evolution of CMB anisotropics that underly the 
calculations which follow. It may be skimmed upon first reading. 

The anisotropy in the CMB is described by small perturbations of the photon distribution function 
around a homogeneous and isotropic black-body. The Boltzmann equation describes the evolution of the 
distribution function /, through Compton scattering with electrons df /drj {r],:x.{r]),p{T])) = C[f], where 
the collision term is written schematically as C[f]. Here rj is the conformal time and p the photon 
momentum. In the absence of spectral distortions, the magnitude of the momentum can be integrated over 
leaving only its directional depenence 7 and the effect of gravitational redshifts on the photon temperature 
perturbation 0. Due to azimuthal symmetry and the decoupling of modes in linear theory, it is convenient 
to decompose the fluctuation in a Fourier or normal mode k into angular moments, e.g. in flat space 
0(77, k, 7) = ^^(— i)^0^Pf (k • 7) with an appropriate generalization to curved spaces ( [Wilson 1983| , |White| 
|fc Scott 1996| ). Here 7i are the direction cosines of the photon momenta. The Boltzmann equation then 
becomes an infinite hierarchy of coupled ordinary differential equations (see e.g. Ma fc Bertschinger 1995 ) 
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where K£ — [1— {P — l)if/fc^]^/^ modifies the angular hierarchy for geodesic deviation in spaces of constant 
comoving curvature K = — 7Jg(l — ilo — ^a) with a Hubble constant of Hq = IQOh km s~^ Mpc""'^. The 
metric perturbations are represented by $ the fluctuation of the spatial curvature in Newtonian gauge and 
^ the Newtonian potential. The collision terms from C[f] are proportional to f = necrra the differential 
optical depth to Compton scattering, where He is the free electron density and ctt is the Thomson cross 
section. 



Scattering by electrons with velocity Vt generates a Doppler effect on the photons. Scattering of 
anisotropic radiation creates a polarization, described by the temperature perturbation in the Stokes 



parameter Q, and governed by a separate hierarchy (Bond & Efstathiou 1984), 
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To complete these equations, we need the baryon Euler equation, which determines the evolution of the 
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baryon velocity, 



vt = --Vb + k^ + T{ei-Vb)/R. (3) 



Finally, the observable anisotropy spectrum follows by integrating over the fc-modes, 

2e+l 1 fdk fc3|9,(7?o,fc)p 
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The interpretation of these equations is quite straightforward. The metric fluctuations feed power into 
hierarchy through the gravitational redshift effects of density dilution ($ in ^ = 0) and potential infall (fc^* 
in ^ = 1). If the optical depth across a wavelength f/fc <C 1, this power flows to higher £ much like a wave 
pulse flows along a string, being concentrated in mode i when ki] ~ £. The critical epoch for this process 
is horizon crossing krj ^ 1 after which ^ > 1 modes can be populated. When the free electron density 
is not negligible then the Compton scattering terms (f terms) become important. Modes with £ > 2 are 
exponentially damped, sealing off the hierarchy and providing a barrier off which the wave pulse reflects. 
The monopole term is not damped at all and the dipolc term is driven toward vt so that the distribution is 
isotropic in the electron rest frame. 

Thus before recombination, f/fc ^ 1, and the photon distribution possesses only the £ = (density) and 
£ = 1 (velocity) modes, which represent a fluid that oscillates acoustically due to photon pressure (see § ^.2[ ). 
Only for very high fc will power leak into the higher £ modes, where it will be exponentially damped. This 
is responsible for the damping tail at small angular scales. An increase in f at late times due to reionization 
also possesses a characteristic signature. For scales inside the horizon at reionization, the power has already 
propagated to high £ where it suffers exponential damping; for larger scales no such damping occurs. Thus 
reionization damps small-scale anisotropics while preserving large-scale anisotropics. We shall discuss these 
behaviors more quantitatively in the next section. 



2.2. Analytic Estimates 

Before turning to the numerical calibration of effects in the damping tail, it is useful to describe them 
analytically to see how they enter in and affect the Boltzmann evolution given above. The two main damping 
processes at work in the CMB are photon diffusion before recombination and rescattering during an epoch 
of late reionization. 



2.2.1. Diffusion Damping before Recombination 

For wavelengths much larger than the mean free path to Compton scattering (fc/f <C 1), the Boltzmann 
hierarchy of Eq. (|^) can be described by the relativistic fluid dynamics of a combined photon-baryon fluid. 
Rapid scattering insures that any anisotropy of the photons in the electron rest frame is vanishingly small 
so that the hierarchy can be truncated at £ = 1 with Qi « Uf,. Even in this tightly- coupled regime, the 
random walk of the photons through the electrons eventually mixes photons across a wavelength of the 
fluctuation ( Bilk 196^ ). Thus, we expect temperature perturbations to be destroyed by this diffusive process 



before the mean free path grows long enough to invalidate the central approximation. This statement is 
only approximately true during the recombination epoch when the mean free path grows so rapidly that it 
approaches the horizon scale and coincides with the diffusion scale. For this reason, we calibrate the process 



numerically in §2.3 
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Formally, we can approximate these effects by expanding the Boltzmann temperature and polarization 
equations ([^) a nd (||) in powers of fc/f ( Peebles fc Yu 197C| ). To lowest order, one obtains the oscillator 
equation ( HSa ) 

^(1 + i?)eo + ^eo = -^(1 + R)^ - + R)^- (s) 

Gravity drives the oscillator by potential infall into and density dilution as the curvature fluctuation <!> 
changes. The baryon inertia in the fluid is described by the relative baryon-photon momentum density ratio 
R and increases the effective mass of the oscillator. Together these effects imply oscillations at the sound 
speed 

cs = , ^ (6) 

^/W+R) 



around a zero point displaced by gravity to (1 + i?)^' for slowly-varying <i> (HSa, HSc, Hu & White f996b) 



To treat the effects of diffusion, one must include higher order terms. An examination of the £ = 1 
photon Euler equation (|^) shows that there are two diffusive effects: viscous damping from the quadrupole O2 
and heat conduction from the relative photon-baryon velocity Qi — Vh ( Weinberg 1972[) . From the expansion 
of the polarization hierarchy, Q2 = Qo = and the quadrupole evolution equation (|^) with Qi 3> 63, we 
obtain the tight-coupling prediction for the quadrupole. 
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Heat conduction may be described by expanding the baryon Euler equation (^ to second order. Let us 
assume a solution of the form Oi cx expi J to drj and ignore variations on the expansion time scale a/ a in 
comparison with those at the oscillation frequency uj. We return to consider this approximation in §3.1. The 
heat conduction equation becomes 
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allowing us to rewrite the photon Euler equation (|l|) as 

iuj{l + i?)ei = fc[eo + (1 + i?)^'] - u^T-^R^Qi 
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The presence of (1 -|- i?)^' again reflects the gravitational zero-point displacement of the oscillator. It is thus 
appropriate to try a solution of the form Oq + (1 + R)'^ oc expi J ujdr]. With this assumption, one obtains 
the dispersion relation for acoustic oscillations 
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From the form of the solution expi J tudr], this gives the damping scale kjj ( Kaiser 1983| ) 
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(11) 



by which acoustic oscillations are damped exponentially as exp[— (fc/Zci))^]. Notice that the diffusion length 
is roughly the geometric mean of the mean free path and horizon length 77 as one would expect of a 
random walk ko ^ ^^r/r/. 
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2.2.2. Diffusion Damping during Recombination 

As the universe recombines, the mean free path and hence the diffusion length of the photons increases. 
As long as the diffusion length is much greater than the mean free path, damping can be described by the 
tight-coupling approximation of the previous section. This is because the mean free path only surpasses 
the wavelength after diffusion has already destroyed the perturbations, resulting in no contributions outside 
of the tight-coupling regime. The approximation thus remains approximately true until quite near the end 
of recombination when the mean free path becomes comparable to the horizon and so the diffusion length 



in §2.3 



(HSa, HSc). This fact explains the reasonable level of agreement between the numerical results we present 



and the tight-coupling approximation. 



The remaining subtlety is that due to the finite duration of recombination, last scattering takes place 
at a slightly different epoch, with a slightly different diffusion length, for each photon. The net effect has 



been approximated (HSa) by weighting the damping by the visibility function re the probability of last 



scattering within drj of 77, 

V{k) = r (re-^) exp{-[k/kD{v)?}. (12) 
Jo 

This "smearing" of the last scattering surface, and the evolution of tend to soften the damping, meaning 
it is not quite the simple exponential one would naively predict. It is however often convenient to define the 
last scattering epoch as t(77,) — 1. 

Notice that the net result depends only on the cosmological parameters of the background. The effect 
of floh^ is simple. Increasing ftoh^ decreases the horizon at last scattering thus monotonically decreasing 
the diffusion length. The dependence on flbh^ is more complicated. Increasing ili,h^ 



1. decreases the mean free path 

2. delays recombination 

3. shortens its duration 

4. speeds diffusion scale growth at recombination 

The first effect tends to decrease the damping length and dominates for low fli,h^. The second effect extends 
the amount of time the photons can diffuse and hence increases the damping length; it dominates at high 
Qbh?. In the limit of instantaneous recombination, the damping function T>(k) attains its sharpest form of 
exp{— [fc/fc£)(r/,)]^}. For the realistic case where recombination takes place over an extended period I?(fc) 
becomes less steep. Both the width of the visibility function and the evolution of kn through it affect this 
drop. Again the baryon dependence of these two effects are in opposition leading to a steepening of the slope 
at both the high and low Q.i,h^ limits. 



2.2.3. Free Streaming 

After recombination photons enter the free-streaming regime. The observer views a temperature 
fluctuation at wavenumber k on the last scattering surface as an anisotropy at multipole moment € « 
fkrs[r]^:), where the constant of proportionality / « 0.98 from numerical fitting (see §2^) and the comoving 
angular-size distance to the epoch ry is 

r,(,7)-|i^ri/2sinh[|i^|i/2(^„_^)], (13) 
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for K < 0. For positively curved universes, replace sinh with sin. The fact that angular size depends 



sensitively on the curvature allows its precise measurement from acoustic features (Doroshkevich, Zel'dovich 



fc Sunyaev 1978, Sugiyama fc Gouda 1992, Kamionkowski, Spergel fc Sugiyama 1994, Hu fc White 1996b) 



These effects are represented formally in the Boltzmann equation (yj) by the transfer of power down the 
^-hierarchy with distance rj — from the last scattering event and the geodesic deviation factors Kg. 
Note that the latter becomes important when the distance is long enough such that the subtended angle 
9 ^ — \/K/k, i.e. smaller than that of a wavelength at the curvature distance. 

The angular-size distance relation may be used to map fc-space inhomogeneities onto ^-space anisotropics. 
For example, the damping function in multipole space is 

V,^V{k = l/fre{r],)). (14) 

There are instances when this mapping fails to accurately describe the streaming process. The projection of 
A:-mode inhomogeneities onto ^-mode anisotropics depends on the viewing angle and is thus not one-to-one. 
In particular, it can take power to larger angles for wavelengths that happen to be viewed with wavevector 
parallel to the line of sight. In this case, the angular separation between the intersections of the flat wavefront 
with the spherical shell at 77* is much larger than {krg)^^. Sharp features in fc-space will thus be blurred in 
£-space and excess power at small physical scales can be aliased into large angular scales. Formally, this is 
reflected by the decomposition of the /c-mode on the sphere and the fact that the solution to the sourceless 
Boltzmann or Liouville equation is just its radial component, a spherical Bessel function in flat space (see 



e.g. Bond fc Efstathiou 1987 ). For the cases we consider, where the A:-space features are broad with no strong 



deviations from scale invariance, the simple approximation of Eq. (0) sufRces 



2.2.4- Reionization Damping 



From the null detection of the Gunn-Peterson effect (|Gunn fc Peterson 1965| ) , we know that the universe 
was reionized at least as early as redshift z « 5. Reionization recouples the photons to the electron-baryon 
plasma. The same process that is responsible for diffusion damping acts to destroy anisotropics during this 
epoch as well. 

During the free-streaming epoch, the effective "diffusion length" is simply the horizon scale. Photon 
trajectories from different temperature regions on the last scattering surface intersect forming the anisotropy 
that is represented by the £ > 2 photon modes. When the universe reionizes, the photons which rescatter 
lose their anisotropy. Note that the isotropic temperature fluctuation that exists above the horizon, where 
trajectories have not yet crossed, does not damp by rescattering. This is reflected in the lack of a f coupling 
term in the £ = mode of the Boltzmann equation (|l|). The £ — 1 mode damps in such a way as to drive 
01 toward Vh so that the distribution is isotropic in the electron rest frame. 

Even in a reionized universe, photons eventually last scatter as the electron density drops due to the 
expansion and the mean free path to scattering exceeds the horizon length. Thus only the fraction e^'^ of the 
photons that did not rescatter contribute to the anisotropy below the horizon at last scattering rj^. Above 
this scale all photons contribute. Thus the rough form of the reionization damping function becomes 

where again the effect of the finite duration of last scattering on rjr can be accounted for by the visibility 
function (see |2.4|). 
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Fig. 1. — Diffusion damping calibration. In the absence of both diffusion damping and gravitational sources, 
the rms temperature fluctuation at recombination (short-dashed line) exhibits simple acoustic oscillations. 
These are mapped onto anisotropies on the sky in a near one to one fashion (solid line). The inclusion of 
diffusion terms in the Boltzmann equation allows for a simple numerical calibration of its effects. 



2.3. Numerical Calibration 

The expressions of the previous section are only approximations, though useful ones. We now turn to 



numerical calibration by solving the Boltzmann equations of §2.1 



Extracting the damping effects from realistic models of structure formation is complicated due to the 
manner in which gravity generates perturbations through the metric fluctuations $ and Vl/ in the model. 
Since the effects discussed above are essentially model- independent, we choose instead to calculate a toy- 
model in which no gravitational effects, beyond the background expansion, are included. Specifically, we 
solve the Boltzmann equations for the photons and baryons with $ = = This includes neglecting the 
self-gravity of the photon-baryon fluid. 

Before recombination, we are left with a pure acoustic oscillation whose behavior is completely 
determined by the initial conditions. For simplicity, we take them to be adiabatic and scale invariant. 



The evolution equations of §2.1 are then solved in the usual way (see e.g. Bond & Efstathiou 1984, Ma & 



Bcrtschinger 1995, Seljak & Zaldarriaga 1996) through recombination to the present. This properly includes 
the effects of diffusion through the last scattering surface and the projection of the fluctuations at last 
scattering onto the sky today. We show an example in Fig. |l| (long-dashed line). 

To extract the diffusion damping behavior, we compare this to a calculation of the same model with 
diffusion damping "turned off" . Speciflcally, we solve the tightly-coupled photon-baryon equation (jsj) up to 
the point where the optical depth to the present, (ignoring reionization) becomes unity. We then free-stream 
the photons to the present by solution of the sourceless Boltzmann equation ('J = $ = f = in Eqs. |l]and^) 
to determine the anisotropy (see Fig. |^, solid line). The ratio of the angular power spectrum of the damped 
to undamped calculation gives the form of the damping function (see also Fig. |^). This also serves to 
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Fig. 2. — Reionization damping calibration. By removing the relative Doppler effect from a reionized 
Boltzmann calculation and comparing the result to the same model (here standard CDM Q,q = 1, h = 0.5, 
Vlbh? = 0.0125) with no reionization, the effects of rescattering damping are isolated. The reionization 
damping envelope is fit by two parameters, the optical depth during reionization and the horizon scale at 
last scattering (see Eq. 24). 

calibrate the angular-size distance relation of Eq. ( |l3| ) through comparison with the mean squared fluctuations 
in Fourier space |0op + |6ip/3 at optical depth unity (see Fig. |l|, short-dashed Hue). By aligning the peaks, 
one extracts the proportionality factor / « 0.98. As discussed in free streaming smears features in the 
fc-space rms spectrum somewhat which explains the slightly smoother actual anisotropy. 

To extract the reionization damping behavior, we compare the no-reionization case to one with some 
arbitrary reionization history. In order to isolate damping effects from the Doppler effect due to the relative 
motion of the baryons with respect to the CMB, we set vi, = Qi during this epoch. For simplicity, we often 
parameterize the reionization as instantaneous at some epoch to some constant fractional level of hydrogen 
reionization Xh though none of our results depend on this simplification. The ratio of the two power spectra 
gives TZj. We show examples in Fig. |^ (solid lines). 



2.4. Fitting Formulae 



It is convenient to fit the numerical calculations of §2.3 for the diffusion damping and reionization 
damping envelopes. Aside from providing a compact summary of the results, this exposes the sensitivity of 
the spectrum to cosmological parameters which will be useful in §pl 
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2.4- 1- Diffusion Damping Envelope 



Since the effect of diffusion damping is determined solely through the Compton mean free path and 
horizon scale, it is dependent on very few cosmological parameters. The Compton mean free path of a 
photon is governed by the baryon density pb cx flbh^- If the present energy density in the radiation is 
fixed, then the horizon only depends on the matter content Q,oh^ before curvature and cosmological constant 
contributions become significant. We assume here that the radiation energy density is fixed by the observed 



CMB temperature Tj ~ 2.728K (Fixsen ct al. 1996) and there exist three families of massless neutrinos 



with = {i/liy/^T^ (we ignore the small correction of pDodelson fc Turner 1992| ) . Thus aside from the 



projection effects from rg, which are sensitive to the curvature and cosmological constant, the damping 
behavior depends only on HqH^ and ilbh^. We have computed T>g as described in §2.3 for 150 models in the 
range 0.02 < Qoh^ < 0.75 and 0.005 < flbh'^ < 0.75 . 

From the tight-coupling expansion, we expect the damping tail to scale as exp[— (^/^c)^]. Including 



the effect of a finite last scattering surface and the conversion from k to £ makes the damping function less 
steep. We find that through the first two decades of damping in power, the function Vg calculated in the 
last section can be approximated as 

= exp[-(£/£cr]. (16) 

The quantities kjjf — In/fe and m are functions of floh^ and flbh^ which are power laws at the extreme ends 
of parameter space (see Fig. ^). Recall that / w 0.98 is obtained by numerical calibration of the projection 



relation (see §2.3). We chose first to fit the flth dependence. A simple two-power law fit 



m = a3(f7fc/l2)'^4[l + (f^^/j2)1.80]l/5 



2^1. 801-1/5 



(17) 



provides a good description of the numerically determined behavior. The coefhcients of this fit are then 
functions of ^^h? . We find that they fit single or double power-laws forms, 



ai = 0.0396(r!o/i')~°'^*[l + 13.6(r!o/i2)"-638]^ 

02 = 1480(f7o/i')~°°''°'^[l + 10.2(r!o/i')°''']~\ 

03 = 1.03(^!o/i')"■°'^^ 

04 = -0.0473(r!o/^2) -0.0639, 



(18) 



Together these fitting functions work to the percent level for 0.02 < fJo/i^ < 0.75 and 0.005 < Jlb/i^ < 0.75 



and improve upon the approximate results of §2.2 and HSc (Eq. E4) 



To complete the description of the damping, we need to express explicitly the conversion of physical 
to angular space variables through the angular-size distance re- The missing ingredient of Eq. ( |l3|) is the 
comoving distance to the last scattering surface rjo — rj^. The horizon scale today can be expressed as an 
integral over the Hubble parameter 770 — J^{a^H)^^da. For Q\ = 0, this has the exact solutions 



1 



Vo 



1 



In 



2-r!o + 2(l-no)'/'(l + acqf^o) 



r!o + 2(l-f^o)i/2(acqr!o)i/' 



tan 



1/2 



2a<i^'(l + Oeq)i/2 



ffo(f^o- 1)1/2 
while for K — the following form. 



tan 



2(1 + acq)(f^o- 1)1/2 



if < 0, 



K > 0, 



(19) 



a/2 



,1/2 



](1 - 0.0841 In f7o). 



(20) 
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Fig. 3. — Diffusion scale calibration. Analytic estimates of koi^oh^ ,^bh^) based on the tight-coupling 
approximation trace the results to reasonable accuracy and explains their general behavior (HSc, Eq. E4). 
The fitting function of Eq. (17) tracks the numerical calibration to better than the 1% level. 



fits the integral over the region 0.1 ^ r2o < 1 and 0.3 < /i to better than 1% accuracy. Here 

acq = 4.17 X lO-%noh^)-\T^/2.728K)^ (21) 

is the scale factor at matter-radiation equality. Finally, the horizon at last scattering where T{ri^) — 1 takes 
the form 

rj, = 2{noH^)-'^^[ia, + a,^)'^^ - Oeq], (22) 



where (HSc, Eq. El) 



-1 = mA8[l + O.Q0124{nbh'^)-°'^^>^][l + biinoh^f^], 
h = 0.0783(r!fc/i2)-o-238[i + 39.5(r!b/i2)0-763]-i^ (23) 

52 = 0.560[l-h21.1(J76/i2)i-8i]-i, 



is a fit to the rcdshift of recombination. 



2.4-2. Reionization Damping Envelope 

Reionization damping depends on two parameters, the total optical depth r and the angular scale 
subtended by the horizon at last scattering during the reionization epoch 9r ^ l^^ . The asymptotic values 
given in Eq. ( |l5| ) are highly accurate and thus we need only search for an interpolating function around l^- 
The following form fits the behavior in TZj to better than 1% for late reionization 



T^i = ' T + (24) 

1 + cix + C2X^ -I- cax-' -I- 042;* 
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with X = £/{£r + 1) and ci = —0.276, C2 = 0.581, C3 — —0.172, C4 = 0.0312. Even the more extreme case of 
early reionization to a low ionization level is described well at the couple of percent level (see Fig. ^). High 
precision in the large optical depth limit is unnecessary since secondary anisotropics dominate in this limit. 

The parameter £r = r0{rir)/rir involves the visibility- weighted horizon at reionization 

^ /d77 7?(fe-n ^ ^25) 
J drj (fe""^) ' 

where the optical depth functions can be obtained by noting that 

T = necrra = 2.304 x 10"^Mpc"i(l - Yp)Q.bh'^{l + zfxh- (26) 

Here Yp « 0.23 is the primordial helium mass fraction and recall that Xh is the hydrogen ionization fraction 
and we assume that helium is not ionized. It is useful to note that for low redshifts Zr <^ 100 and constant 
ionization fraction, the optical depth may be integrated analytically to give 



T = 4.61 X 10-2(1 _ Yp)xJ^ {2 - 3f^o + (1 + n^ZrY'^i^oZr + 317o - 2)| (27) 
when A = and 

T - 4.61 X 10-2(1 - Yp)xJ^ {[1 - r!o + W + ZrfY'^ - 1} (28) 



h 

when K = Q. For higher redshifts, the contribution of the radiation to the expansion rate can make a few 
percent or greater correction. 



3. Cosmological Information 



Armed with the calibration of the effects of diffusion and reionization damping, we can now examine 
the information, both on cosmological parameters and models for structure formation, embedded in the 
small-scale anisotropy spectrum. We begin with a discussion of the assumptions that render diffusion and 



reionization damping model-independent for most models of structure formation (§3.1). By removing the 
effects of damping in such models, one uncovers striking signatures that contain essential information on the 



nature of fluctuations in the early universe (see also |Hu fc White 1996bD . For illustrative purposes, we often 
employ variants of the standard CDM model, scale invariant initial adiabatic fluctuations A:'^|^'(0, k)\'^ =const. 
in an r^o = 1, /i = 0.5, fib/i^ — 0.0125 universe. 



Baryon drag, which enhances alternate acoustic peaks (§3.2), can help separate adiabatic from 
isocurvature fluctuations, an important step in distinguishing inflationary models from cosmological defect 



models (see |Hu fc White 1996q , |Hu, Spergel, fc White 1996]) . It also probes the gravitational potential at 
last scattering. The underlying amplitude of the oscillations extracts information on the evolution of the 



gravitational potentials at horizon crossing through the potential envelope (§3.3) 



The diffusion damping and reionization damping envelopes are themselves interesting because they 
provide essentially model-independent information about cosmological parameters, mainly the curvature of 



the universe (§3.4) and the epoch and extent of reionization (§3.5). In this section, we systematically treat 
these applications of the results from the damping calibration in §^. 
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3.1. Model Assumptions 



We begin by examining the conditions under which the diffusion damping and reionization damping 
enveiopes are modei-independent to expose general guidelines for their use. 

Only acoustic oscillations are damped by diffusion. This leaves untouched e.g. offsets in the zero point of 
the oscillations or anisotropies generated between the last scattering surface and the observer. In the former 
case, the — ^' offset provided by the potential is not damped because it represents gravitational redshifts 
which are picked up by the photons even as they diffuse in and out of potential wells. The baryons provide 
an inertia to the photon-baryon fluid which further offsets the oscillation. The Compton drag of the baryons 
on the photons increases the photon temperature inside gravitational potential wells by —R'^ leading to a 
zero-point shift that is not damped by diffusion for similar reasons. Together the redshift and drag effects 
explain why in the estimates of §2.2, it is Oq + (1 + i?)^ which suffers damping and not Oq- 



The time evolution of the potentials causes a shift of order ^/k^ [see Hu fc White 1996b , Eq. (25)]. If 
i? ^ |'I'//c^\E'|, it is negligible in comparison to baryon drag. Generally ^E" varies on the order of an expansion 
time such that |5'/fc^\E'| = 0[{krj)~'^] ^ 1 for scales well inside the horizon at last scattering: krj ^ 1. Mixed 
terms of the order R"^ /k also exist but are again generally smaller than the R^ term. Since the intrinsic 
acoustic amplitude is of order the gravitational potential at sound horizon crossing \E'(A;,r~^), the diffusion 
damping signature dominates over the undamped term if 



and 



> R 



>1, 



(29) 



(30) 



*(fcD,?7*) 

which are generally satisfied by models whose potentials do not grow significantly well within the sound 
horizon. Notice that no assumption of coherence in the oscillation is necessary (Magueijo et al. 1996| ). 



In principle, there is also a model-dependent effect since in the discussion above we have implicitly 
assumed a two-step process: first the acoustic oscillations are formed and then they are damped. This is 



generally called in the literature a "passive" approximation ( Albrecht et al. 1995 ). If the model possesses a 
strongly time- varying potential inside the horizon, the underlying acoustic oscillations could still be forming 
as the diffusion length overtakes the wavelength. Usually this is a small effect, since most of the damping 
occurs at the instant of recombination so that the fluctuations generated during this short time are small. 

Finally anisotropies generated between recombination and the present could be larger than the intrinsic 



acoustic signal, especially in the damping tail. This could occur due to the linear (Kaiser 1984) and non- linear 



(Sunyaev & Zel'dovich 197C, Vishniac 1987) Doppler effects in a reionized universe or time variations in the 
potential along the line of sight (Sachs & Wolfe 1967, Rees & Sciama 1968, Kaiser & Stebbins 1984). In models 



such as cold dark matter (CDM) with a near scale-invariant spectrum of adiabatic initial fluctuations, this is 
not a worry. The lack of excessive small scale power in the model makes early reionization and/or small-scale 
non-linearities that are responsible for such effects unlikely ( |Ceballos fc Barcons 1994 , $eljak 1996| , |Hu fc| 



White 1996a| ). 



These types of considerations also apply to the reionization damping function TZi calculated in § ^.3| . 
By construction, this function isolates the rescattering damping effect during reionization and ignores any 
secondary effects that may regenerate fluctuations. Again, the Doppler effect due to the relative velocity of 
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Fig. 4. — Baryon drag and its potential dependence. Baryon inertia in the fluid displaces the zero point of 
the temperature oscillations leading to alternating peak heights as a function of scale at last scattering. The 
magnitude of the displacement is i?*^(77*), and by removing it the monotonic variation of heights due to the 
potential envelope is uncovered (upper panel). The fractional effect is of order R^,^{r|^,, k)/'i'{0, k) and can 
be adequately described by the matter transfer function T{k) (lower panel). The model here is CDM with 
nQ = l,h = l and Vtbh? = 0.025. 



the electrons with respect to the CMB can regenerate fluctuations significantly if both the peculiar velocities 



and the optical depth are large. We examine this effect more closely in §3.5 



In summary, the damping function I?£ accurately describes the model-independent damping of acoustic 
oscillations and the reionization damping function TZg does the same for the rescattering damping oi primary 
anisotropies. In models such as CDM with no excess small scale power and hence relatively late reionization 
and small secondary effects, their behavior will be clearly manifest in the observable spectrum. In models 
where this is not true, it merely describes the behavior of a component of the total anisotropy and other 
effects must be taken into account to extract the information embedded in the observed anisotropy. 



3.2. Uncovering the Baryon Signature 



Baryons create a distinct acoustic signature due to the drag effect discussed in §3.1. By providing inertia 
to the fluid, they enhance compressions over rarefactions inside potential wells. Aside from providing a means 
to measure the baryon content, it also distinguishes between between the two phases through the difference 
in peak amplitudes between successive peaks. In turn this distinction provides one of the most striking and 
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robust ways to distinguish adiabatic inflationary fluctuations from their isocurvature counterparts, generated 



perhaps by cosmological defects ( Hu fc White 1996b , Hu, Spergel, fc White 1996 ). Unfortunately damping 



and projection effects serve to obscure this signal. By deconvolving these effects with the results and methods 



of §2.3-2.4, one can uncover this important signature. 



Let us first examine the intrinsic effect. In Fig. ^, we show an example from a solution of the tight- 
coupling oscillator equation Eq. ^ under the metric fluctuations of a,n flo — 1, h — 1, flth^ — 0.025 
CDM model. Displayed is the effective temperature fluctuation of the peaks (triangles), connected by the 
full function to guide the eye. To demonstrate that the alternating height effect is due to baryon drag, 
we add i?,^'(?7*,fc) to each peak (squares), where i?, = Riv*)- Notice that this eliminates the alternation, 
leaving the peak heights to smoothly vary in a manner described by the "potential envelope" discussed in the 
next section. Since the intrinsic amplitude of the oscillations is of order the potential before sound horizon 
crossing, the fractional effect is of order i?*^'(ry, , fc)/4'(0, fc). 

Since the effect depends on the potential at last scattering '^{r]^,,k), it also provides a probe of the 
matter fluctuations at that epoch. Under the CDM scenario, the potential does not evolve significantly 
between recombination and the present so that the baryon drag effect also reflects the matter fluctuations 
today. The fractional effect becomes i?**(r7*, fc)/\I'(0, fc) « R^'^{r]a, fc)/*(0, fc) = R*T{k), where T(fc) is the 



matter transfer function ( Bardeen et al. 1987 ) 



nq) = ^^^^^^^^^[1 + 3.89g + {le.lqf + (5.46g)3 + i6.71q)Y'^\ (31) 

with q = (fc/Mpc"^)(r^/2.7K)2(Oo/i^)"^ In Fig. | (lower panel), we show that RT{k) accurately tracks 
the effect and provides a potential consistency check with large scale structure today. As we shall see in the 
next section, the fall of the fractional baryon drag effect and the rise of the potential envelope are intimately 
related through the matter-radiation equality epoch. 

The magnitude of the baryon drag effect in the observable anisotropy spectrum is reduced by inclusion 
of the dipole term and smoothing by projection, but mainly by diffusion damping. If the baryon content 
is low, the intrinsic magnitude of the effect is small and diffusion damping may cause the peak heights to 
monotonically decrease rather than alternate (see Fig. Given the calibration of the diffusion damping 
behavior in §2.3, we can invert this flltcr. In Fig. ^, we show that multiplying the spectrum by 'D'^^ uncovers 



the alternating peaks even for fl^h^ signiflcantly lower than the standard big bang nucleosynthesis prediction. 
In practice, removing the damping behavior will require knowledge of flgh^ and fl(,h^ either from external 



measurements or consistency checks (see Hu & White 1996b) as well as measurement of the curvature from 
the CMB. 



3.3. Determining the Potential Envelope 



Gravitational potential perturbations drive acoustic oscillations affecting their amplitude and phase. The 
effect on the phase can be used to uncover information about the origin of fluctuations in an inflationary 



epoch or phase transition (HSb, Crittenden & Turok 1995, Hu & White 19961). Here we treat their effects 



on the amplitude of the intrinsic oscillations, unobscured by the presence of diffusion damping. This can be 
obtained from an observed spectrum by the techniques of §0 and is also useful for constraining the curvature 



(see gSJ) 

As an example of the driving process, let us consider the case of adiabatic fluctuations. The self- 
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Fig. 5. — Uncovering Baryon Drag in a low baryon universe. Diffusion damping obscures the baryon drag 
signal especially in a low baryon universe (here n^h^ ~ 0.075 in an otherwise standard CDM model). 
Employing the numerical calibration of the damping tail, we recover the alternations. 



gravity of the photon-baryon fluid drives its own oscillations through a feedback mechanism at sound horizon 
crossing. Photon pressure prevents gravitational collapse inside the sound horizon leading to a decay in the 
self-generated gravitational potential. This decay is timed such that it leaves the oscillator in a highly 



compressed state leading to correspondingly large amplitude acoustic oscillations (see Hu & White 1996b 
for further description). If the self-gravity of the photons and baryons dominate, the amplitude of the 
oscillation is enhanced from gravitational redshifts by — 2^* which combined with the Sachs- Wolfe effect 



(Sachs & Wolfe 1967) of yields a net result of — l^t, i.e. the amplitude of the oscillation should be 5 times 
the large-angle Sachs- Wolfe plateau. Inclusion of neutrinos and the matter-radiation transition modify this 



result to 5(1 -I- 15/1/)"^, where the neutrino density fraction is — p^/iPv + P-y) ( HSc , Eq. B9). This driving 
effect only operates if the self-gravity of the photon-baryon fluid dominates at sound horizon crossing. Large 
scales cross the sound horizon in the matter-dominated epoch and do not suffer this effect. Thus the scale 
that crosses the horizon at matter-radiation equality /ccq marks the transition between the two asymptotic 
regimes. 

The critical scale /ccq = (2iloi?o/'^oq)^^^ oc Vl^h?' provides the CMB with sensitivity to the parameter 
Q,Qh? . This is similar to the more familiar effect of equality on the matter power spectrum [see Eq. (31)] but 
note that fluctuations increase rather than decrease upon crossing fceq. Fig. ^ shows that the potential 
envelope that governs the amplitude is indeed a function of k/k^q oc k/floh^. Potentially, this effect 
can also probe the neutrino mass through its effect on 
Bertschinger 1995| , iDodelson, Gates fc Stebbins 1995| ). 



(see Seljak &: Bertschinger 1993, Ma fc 



The remaining subtlety is that the presence of baryons makes acoustic oscillations decay adiabatically. 
Notice that the tight-coupling equation (||) describes an acoustic oscillator with effective mass of (1 + i?). The 
adiabatic invariant for such an oscillator is the energy/frequency. This requires that temperature fluctuations 
decay as (1 -I- R)~^^^ and dipole or Doppler contributions to decay as (1 -I- R)^^/^. The amplitude of the 
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Fig. 6. — Potential Envelope. Decay of the potential due to the self gravity of the photon-baryon fluid drives 
the oscillator. Comparing two CDM models with differing matter to radiation ratios floh^, we see that the 
oscillations are multiplied by an envelope that depends on the equality scale fccq oc floh^. 



potential envelope thus gains a baryon dependence set by the value of R at recombination (HSc). 
The full potential envelope in power can be roughly described by 

Vt. = l + Acxpi-lA£,Ji), (32) 

for a scale-invariant spectrum. Here ^eq = keqre and the amplitude A is fixed by the asymptotic expression 



(33) 



where we have combined the temperature and Doppler effects in quadrature. Tilting the primordial spectrum 



produces an analogous tilt in Ve- The integrated Sachs- Wolfe effect (Sachs & Wolfe 1967) in open and A 
models, also gives Vi large-angle contributions (see also Hu, Sugiyama fc Silk 1996). 



We show an example in Fig. ^. The upper curves show a calculation with the effects of diffusion damping 
removed compared with the potential envelope of Eq. (|3^) . Note that the form of the envelope roughly traces 
power in the fluctuations. The bottom curves show how diffusion damping obscures the signature and tests 



the damping calibration of §2.3 in a realistic context. By multiplying the undamped calculation by one 
regains, to reasonable accuracy, the result of a full CDM calculation encorporating diffusion damping. 

Thus the obscuring effects of diffusion damping can be removed to extract the potential envelope of 
acoustic oscillations. This provides information on the evolution of the metric fluctuations as they cross the 
sound horizon which may help unravel information about the nature of such fluctuations in the general case 
and the scale of matter-radiation equality in the adiabatic case. 
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Fig. 7. — Uncovering the Potential Envelope. The potential envelope is obscured by diffusion damping. By 
numerically removing the damping, one sees that the intrinsic fluctuations follow the analytic estimates of 
Vi reasonably well. By multiplying by the numerically-calibrated damping function T>J, one recovers the 
form of the full calculation even at very small angles. The model here is standard CDM. 



3.4. Constraining the Curvature 

The angular scale of diffusion damping ijj provides a clear feature by which a classical angular-size 



distance test of the curvature can be made by comparison with the corresponding physical scale kjj (Hu & 



White 1996a). In models with simple acoustic peak features which can also be used for this test, this provides 



a consistency check on the curvature, important if the baryon content or thermal history of the universe is 



unknown or anomalous (Hu & White 1996b). In models where the peak signature is more complicated or 



non-existent (Albrecht et al. 1995), it may serve as the primary means of measuring the curvature. 



In principle, the curvature is constrained by the simple absence or presence of small scale power. In 
an open universe, geodesic deviation moves the diffusion tail in angular space to smaller angles leading to 
more power on small scales. In practice, its application is complicated by secondary effects in the foreground 
and lack of a priori knowledge about the intrinsic amplitude of fluctuations before damping. The former 
is unlikely to be an obstacle in models with no strong non-linearities at small scales, in which the acoustic 
signal from recombination is the dominant contribution to the anisotropy. 

Lack of knowledge of the intrinsic amplitude of oscillations limits the precision by which the curvature can 
be measured from the damping tail. The intrinsic amplitude is given by the potential envelope Vg discussed 
in the last section. Given that diffusion damping is exponential in £, it takes only a rough estimate of Vi to 
yield interesting constraints on the curvature. Furthermore if Vg is a slowly varying function compared with 

measurement of the power on several scales in the damping tail can remove the ambiguity. 

We show an example in Fig. ||. Here we assume that the underlying spectrum is of standard CDM which 
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Fig. 8. — Constraining Qq with the damping tail. By measuring the anisotropy power in at some scale i 
in the damping tail (here averaged over 10% in £) and comparing it to a reference scale (here £ — 2), one 
determines the ratio of intrinsic powers Vi /V2 before damping necessary to reproduce the observation (here 
fio = 1 in standard CDM). Since this is a strong function of the assumed fio, only order of magnitude 
knowledge of the model-dependent intrinsic power is needed (e.g. square, estimated from Eq. (32)) to reject 
values of f2o- Multiple measurements in the damping tail largely removes this ambiguity (curve intersection). 
For simplicity, we have fixed h = 0.5, fli,h^ = 0.0125 and Ha ~ 0. 



sets f^o = 1. By comparing the power at some scale i in the damping tail to some reference scale, here £ — 2, 

£{£ + l)Ci Vj Vi 



6C2 P| V2 



(34) 



and by using the fitting formula for T^i of Eq. (|l7|), one can determine the intrinsic ratio of power VijV^ as a 
function of f^o needed to reproduce the measurement. We have ignored here the suppression of power from 
IZg as it is generally negligible for our purposes here. Because the damping multipole is a strong function 
of r^o, the amount of intrinsic power required increases steeply with fig- Thus even the crude estimate of the 
CDM potential envelope of Eq. (^2|) is more than sufficient to distinguish between interesting values of JIq 
(see Fig. ^, square). As the slope of the curves reflect, the further into the damping tail one can measure, 
the more powerful the test becomes. Of course for l':^ Id, the signal also drops exponentially and hence is 
difficult both to measure and separate from secondary effects. 

By measuring more than one scale in the damping tail, one obtains a consistency check on the curvature 
constraint. If the i dependence of is weak, as is the case for CDM-like scenarios (see Fig. |^), then the 
predictions for the intrinsic power must intersect near the actual value of fio • This implements the damping- 
tail shape test proposed in Hu fc White 1996b| to remove the model dependence of the curvature constraint. 



Even if only upper limits exist from CMB measurements at small scales, lower limits on fio can be obtained 
with reasonable assumptions on the amount of intrinsic power in small scale fluctuations. 
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I 

Fig. 9. — Rcionization damping in standard CDM. Damping described by the envelope TZi is the main effect 
of late rcionization in CDM type models. Hence employing either the numerical calibration of TZ^ and the 
fit to it from Eq. (24) to filter the results of a standard recombination (SR, no rcionization) calculation 
approximate the full calculation to better than 1% in power. The scatter at low £ is a numerical artifact 
from finite sampling of the integral in fc-space (see Eq. (4)). 



3.5. Examining Reionization 

Even late reionization produces potentially observable consequences for precise measurements of the 
CMB. In a standard CDM model, the optical depth ranges from 1-3% between 5 < < 10 leading to a 
2-6% effect in the anisotropy power spectrum. For these low optical depths, it is likely that the main effect 



of reionization is the rescattering damping calculated in §2.3, In this case, two cosmological quantities are 
potentially extractable from the spectrum, the total optical depth and the horizon size at last scattering 
during the reionized epoch. In practice, extracting accurate results will be hampered by cosmic variance at 
large angles and the close degeneracy between changes in the spectrum due to the normalization and late 
reionization at small angles. 



In Fig. |[ we show how well the numerical calibration of §2_^ and the fitting formula of §2^ Eq. ( p^ ) 



reproduce the full effect of late reionization. The accuracy achieved is always better than a percent in 
power and increases toward small scales where the reionization signal is the largest. With the high precision 
achievable by the next generation satellite experiments, it is conceivable that the CMB spectrum can probe 
even such relatively late reionization. 

If reionization occurs earlier, such that the optical depth is higher and/or non- linear effects dominate, 
then its effect on the CMB can be even more significant. Fluctuations are not only erased but also 
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Fig. 10. — Reionization and the Doppler Effect. For early ionization, the Doppler effect due to the relative 
electron-photon velocity can regenerate fluctuations around the horizon scale at the last scattering epoch. By 
comparing the standard recombination (SR) result filtered by reionization damping TZ"^ to the full calculation, 
we can uncover such effects. 



regenerated. As an example, consider the Doppler effect from the relative velocity of the electrons with 
respect to the CMB generated as the baryons fall into dark matter potential wells. Its effect peaks near the 
horizon at last scattering due to competing effects. Velocity flows are only generated inside the horizon. Yet 
on small scales, photons last scatter against many crests and troughs of the velocity perturbation leading to 



a strong cancellation damping of the Doppler effect (Kaiser 1984). By employing the rescattering damping 
function TZi from § 2.3 , we isolate this effect in Fig. |1C|. For the higher optical depth cases, the Doppler effect 
is clearly apparent as an excess of fluctuations from that predicted by TZi. On scales much smaller than the 



horizon at the last scattering epoch, simple analytic approximations exist for this effect ( Kaiser 1984 , Hu fc 
White 1996a). In a CDM model where the optical depth is likely to be r < 1, such small scale effects are 



masked by larger primary anisotropics until well into the damping tail. 

More complicated rescattering damping can occur if the reionization is patchy. Although one cannot 
directly apply the results of our damping calibration to this case, basic elements uncovered such as the 
dependence of damping on the horizon scale can be applied to this case as well. Non-linear effects can also 
create fluctuations through the Doppler effect but these are generally small in a model like CDM without 



excessive small scale power (but see Aghanim et al. 1996) 
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4. Conclusions 

Prospects for measuring the small scale CMB anisotropy spectrum are bright, especially in light of the 
approval of two new satellite missions, MAP from NASA and Cobras/Samba from ESA, and the funding of 
ground based interferometers. If foregrounds, systematic and secondary effects are small or can be removed, 
and the inflationary CDM model is correct, much cosmological information can be extracted from the 



damping tail of CMB anisotropics (see e.g. BcrsancUi et al. 1996 ). Despite the enormous success of this 



model however, it is quite possible that what is found there will come as a surprise to the current orthodoxy 



in cosmological modeling. In preparation for this possibility we have here, and in Hu & White 19961 , 
attempted to construct the spectrum out of fundamental physical effects that are likely to be the elements 
in any future model which successfully explains the observations. 

The basic elements uncovered here represent a series of numerically calibrated transfer functions that 
describe the linear processing of fluctuations: the diffusion damping envelope, the reionization damping 
envelope, the potential envelope and the baryon drag modulation. The anisotropy spectrum is not merely 
a snapshot of conditions on the last scattering surface. Rather it is a dynamic entity that bears the mark 
of fluctuations before horizon crossing through the acoustic phase ( Hu fc White 1996b| ), at horizon crossing 



through the potential envelope, at last scattering through baryon drag, and after last scattering through the 



large-angle potential envelope ( Sachs fc Wolfe 1967 ) as well as the effects of reionization. Within the present 
framework of model possibilities, this view of its structure also creates a system of consistency checks by 
which we can verify model assumptions, such as the inflationary or cosmological defect origin of fluctuations, 
before proceeding to measure cosmological parameters and details of the model. 
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